LBL-35097 






> 

o 
a: 



X 

£ 



COLOR, SPIN AND FLAVOR DIFFUSION IN QUARK-GLUON PLASMAS [j 



H. HEISELBERG 



Nuclear Science Division, MS 70A-3307 

Lawrence Berkeley Laboratory, 

University of California, Berkeley, California 94720 



January, 1993 



'This work was supported in part by DOE grant No. DE-AC03-76SF00098, and the Danish Natural Science Research Council. 



COLOR, SPIN AND FLAVOR DIFFUSION IN QUARK-GLUON PLASMAS 



H. Hciselbcrg 

Nuclear Science Division, MS 70A-3307 

Lawrence Berkeley Laboratory, 

University of California, Berkeley, California 94720 

(Received January 1994) 



In weakly interacting quark-gluon plasmas diffusion of color is found to be much slower than the 
diffusion of spin and flavor because color is easily exchanged by the gluons in the very singular forward 
scattering processes. If the infrared divergence is cut off by a magnetic mass, m mag ~ a s T, the color 
diffusion is D co i or ~ (a s ln(l/a s )T) _1 , a factor a s smaller than spin and flavor diffusion. A similar 
effect is expected in electroweak plasmas above Mw due to W exchanges. The color conductivity 
in quark-gluon plasmas and the electrical conductivity in electroweak plasmas are correspondingly 
small in relativistic heavy ion collisions and the very early universe. 

PACS numbers: 12.38.Mh, 12.38.Bx, 52.25.Dg 



Transport of color degrees of freedom in a quark-gluon 
plasma has recently been found to be infrared sensitive 
0] and thus differ from other transport processes as vis- 
cous and thermal flow, stopping and electrical conduction 
PH as well as energy degradation W. Color does not 
flow easily due to the transfer of color (color-flip) in the 
exchange of a colored gluon in the very singular forward 
collisions. As will be discussed here, this suppression of 
the color flow also applies to the color diffusion and color 
conductivity, which are infrared sensitive like the quark 
and gluon quasiparticle relaxation rates MM- These re- 
sults do not only apply to QCD plasmas consisting of 
quarks, antiquarks and gluons, as existing in the Early 
Universe and searched for in relativistic heavy ion colli- 
sions at Brookhaven (AGS and RHIC) and CERN (SPS 
and LHC). As will be shown below it is a general mecha- 
nism for most non-abelian gauge theories. For example, 
in an electroweak plasma consisting of leptons, photons 
and electroweak bosons at temperatures above M\y ~ 80 
GeV, the exchange of W ± provides a charge exchange 
analogous to the color exchange which results in corre- 
spondingly poor electrical diffusion and conduction. 

Earlier studies of transport processes in relativistic 
quark-gluon and electron-photon plasmas found that the 
effect of Landau damping effectively led to screening of 
transverse interactions and gave the characteristic relax- 
ation rates in transport processes. Transport coefficients 
for weakly interacting electron-photon and quark-gluon 
plasmas for both thermal plasmas &fl] as well as de- 
generate ones |7J] were calculated to leading logarithmic 
order. Generally the transport relaxation rates have the 
following dependence on interaction strength 
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l/r tr ~a 2 ln(l/a s )T. 
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However, the quark and gluon quasiparticle damping 
rates, l/r p , were not sufficiently screened by Landau 
damping for non-vanishing quasiparticle momentum, p, 
and depends on an infrared cut-off, m ma g — a s T, so that 
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to leading logarithmic order. Since the quasiparticle de- 
cay rates are not measurable transport coefficients the 
infrared sensitivity was not considered a serious problem. 
However, it was recently discovered Q that diffusion of 
color in some abstract color space suffered from the same 
infrared divergence which led to the same color relaxation 
rates as the quasiparticle damping rates. 

We will describe the two kinds of transport processes 
by calculating the flavor, spin, and color diffusion coef- 
ficients in a quark-gluon plasma within the Boltzmann 
kinetic equation 
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Here p^ and £j are the quasiparticle momentum and ( 
ergy respectively, rii(pi) the quasiparticle distributi 
function, and F the force on a quasiparticle. The r.h 
is the collision integral for scattering particles from i: 
tial states 1 and 2 to final states 3 and 4, respectivf 
with matrix element squared |Mi2^34| 2 summed over 
nal states and averaged over initial states. The (1 ± 1 
factors correspond physically to the Pauli blocking of 
nal states, in the case of fermions, and to (induced ( 
stimulated emission, in the case of bosons. v 2 is the si 
tistical factor, v 2 — 16 for gluons and v 2 = YINf : 
quarks and antiquarks. For scattering of quarks of c 
ferent flavor 
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quark-gluon and gluon-gluon interactions are just 9 
and (9/4) 2 times stronger respectively near forward sc; 
tering. In a medium this singularity is screened as giv 
by the Dyson equation in which a gluon self-energy n^ 
is added to the propagator 
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(we refer to B for details on separating longitudinal a 
transverse parts of the interaction) where the longitu* 
nal and transverse parts of the self-energy in QED a 
QCD are for u,q <C T given by 
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where x = Lo/qv p and v p — c for the relativistic plasn 
considered here. The Debye screening wavenumber 
thermal QCD is q 2 D = g 2 (2N + N f )T 2 /6 where N = 3 
the number of colors, Nf is the number of quark flavo 
T the plasma temperature and /j, q the quark chemii 
potential. We refer to || for a detailed comparison 
QED plasmas. In the static limit, IIi(w = 0, q) = q 
and the longitudinal interactions are Debye screens 
For the transverse interactions the self-energy obeys t 
transversality condition q^W^ = 0, which insures th 
the magnetic interactions are unscreened in the sta 
limit, ny(w = Q.q) = 0. It has therefore been si 
gested that the transverse interactions are cut off beh 
the "magnetic mass" , m mag ~o 2 T, where infrared div> 
gences appear in the plasma [ |I0| . However, as was she 
in P,B|, dynamical screening due to Landau damping 
fectively screen the transverse interactions off in mt 
transport problems at a length scale of order the Del 
screening length ~ 1/gT as in Debye screening. Nc 
ertheless, there are three important length scales in t 



quark-gluon plasma. For a hot plasma they are, in in- 
creasing size, the interparticle spacing ~ 1/T, the Debye 
screening length ~ 1/gT, and the scale l/m mag ~ l/g 2 T 
where QCD effects come into play. A weakly interacting 
QCD plasma and its screening properties is very similar 
to a QED plasma if one substitutes the fine structure 
constant a s — g 2 /4ir by a = e 2 /Att ~ 1/137, the gluons 
by photons and the quarks by leptons with the associated 
statistical factors || . 

Let us first consider a quark-gluon plasma where the 
particle flavors have been separated spatially, i.e., the 
flavor chemical potential depends on position, /ii(r). In 
a steady state scenario the quark flavors will then be 
flowing with flow velocity, Ui. For simplicity we take the 
standard ansatz for the distribution functions (see, e.g., 
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The expansion is valid near equilibrium where \ii and 
therefore also u^ is small. It gives two terms: the equilib- 
rium distribution function n\ = (exp[(e p — /Zi(r)/T]) ± 
1) _1 and the deviation from that. In general the devi- 
ation from equilibrium has to be found self-consistently 
by solving the Boltzmann equation. However, as in the 
case of the viscosity || , we expect the ansatz (||) to be 
good within few percent to leading logarithmic order. 
The flavor diffusion coefficient, Dfi avori defined by: 
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is given in terms of the flavor current j^ and the gradient 
of the number density Pl = £ n°(p;) = ^r 3 3£(3)/47r 2 
of a particular flavor i. From (pj) we find 
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The density gradient, Vp; = ^^iJ2 p (dn^/ds p ) can be 
found by solving the Boltzmann equation. Linearizing in 
Ui we obtain 
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It is most convenient to choose the plasma center-of-mass 
system where one flavor is flowing with velocity u 4 and 
the others with velocity u 2 = — Ui/(7V/ — 1). The num- 
ber of scatterers is then v 2 = 12{Nf — 1). Equivalcntly, 
one can conveniently include the first flavor so that the 
number of scatterers is v 2 = 12JV/ but u 2 = 0. In steady 



state the gluons will not move in the c.m.s., i.e. u 2 = 
for quark-gluon scattering. Since the flavor is unchang 
in the collisions 113 = Ui and 114 = u 2 . 

To leading logarithmic order the singular interacti 
near forward scattering allows us to expand around q 
0, where q = pi — P3 = P4 — p 2 is the momentum trans 
in the collision. Multiplying both sides of ( pd| ) by pi a 
summing the Boltzmann equation reduces to 

niV/*i = -u 1 |i^j Y, n° 1 n 2 (l-n° 3 )(l±n'i) 
q,pi,P2 
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where we have used the antisymmetry of the r.h.s. 
coordinate change pi — > p 3 so that pi — > q/2. The r.b 
collision integral of Eq. (12|) is now straightforward 
evaluate to leading logarithmic order when the screeni 
is properly included (see also Refs. y^Jj). We find 
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where £, is the Rieman zeta-function. The term " 
arises from quark-gluon scatterings and the JV//6 frc 
quark-quark scatterings. This result is similar to the v 
cous, thermal and momentum relaxation rates becai 
the collision term contains the same factors of moms 
turn transfer: the singular q~ 4 factor from the mati 
element squared and the suppressing q 2 factor becai 
the quark flavors lose little momentum in forward sc; 
terings. Including screening, g~ 4 — * (q 2 + II^t)" 
where effectively TVl,t ~ q%, and integrating over n 
mentum transfer, d 2 q, gives the leading logarithmic tei 
HT 2 /q 2 D )c\n(l/a s ). 

Subsequently, let us consider the case where the par 
cle spins have been polarized spatially by some magne 
field J 11 1, i.e., the spin chemical potential depends 
position, /v(r). With the analogous ansatz to (g) : 
the distribution functions with \i a instead of //,, we fi 
the spin current j CT = \i a n a . Linearizing the Boltzma 
equation we find 

Jiv m ^- Vl =2^ 2 Y, <nl{l-nl){l±nl) 
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where M** and M^ are the amplitudes for interacti 
with and without spin-flip respectively. Without spin-f 
the usual factor q as in flavor diffusion appears. Wi 
spin-flip, however, u 3 = u 2 and U4 = U! and the fad 
(Pi — P2) appears. Due to Galilei invariance both ten 
are necessarily proportional to the relative flow, (ui — u 
The transition current can be decomposed into int< 
actions via the charge and the magnetic moment by t 
Gordon decomposition rule 
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where only the latter can lead to spin-flip. We notice that 
the spin-flip amplitude is suppressed by a factor pf — pi 
which leads to a spin-flip amplitude suppressed by a fac- 
tor q 2 . We then find that the spin-flip interactions do not 
contribute to collisions to leading logarithmic order and 
the collision integral is similar to that for flavor diffusion 
evaluated above. Consequently, the corresponding quark 
spin diffuseness parameter becomes 

D^ = D flavor . (16) 

Gluon spin diffusion is slower by a factor 4/9, due to the 
stronger interactions, and by another factor 4/9, due to 
differences between Bose and Fermi distribution function, 

i.e. Di 9) ~ (4/9) 2 ^ 9) . 

Finally, let us, like for the spin diffusion, assume that 
color has been polarized spatially given by a color chem- 
ical potential, /U c (r). The basic difference to flavor and 
spin diffusion is that quarks and gluons can easily flip 
color directions in forward scattering by color exchanges, 
i.e., one does not pay the extra q 2 penalty factor as in 
the case of spin-flip. Consequently, the color-flip inter- 
actions will dominate the collisions since they effectively 
reverse the color currents. The Boltzmann equation thus 
gives us an analogous result to Eq. (JL4|) replacing spin 
by color where the color-flip amplitude now dominates. 
The flow velocity of the scatterers, u,, i=1...4, depends 
on what color combination of the scattering quarks, an- 
tiquarks and gluons. However, in c.m.s. the scatterer 
has vanishing flow velocity, U2 = 0, on average. Likewise 
the final velocities will be zero on average. Multiplying 
both sides with pi and summing the Boltzmann equation 
reduces to (c.f., Eq.(|l4j)) 

niV(ii ]C = — Ui7w 2 /J n\n\{l — n®){l ± n^) 
q,pi.P2 

x |Mi2_>3 4 | 2 (pi - P2) 2 S(e 1 +e 2 -e 3 - e 4 ), (17) 

where we have used the antisymmetry by interchange of 
Pi - * P2- The matrix element entering in ( |17| ) is now 
averaged over all color combinations. 

The transverse interactions actually diverge for small 
momentum and energy transfers even when integrating 
over energy transfers, i.e. dynamical screening is insuf- 
ficient for obtaining a non-zero color diffusion coefficient 
like for the quasiparticle decay rates in QCD and QED 
plasmas (see ||). Concentrating therefore on the lead- 
ing contribution from transverse interactions for small 
jU = uj/q, where n^ — 1(^/4)5^^, we find to leading 
order 
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where we have introduced an infrared cutoff, A. The 1 
per limit on momentum transfers, ~ T, actually con 
from the distribution functions in Eq. ( [17] ) but it dc 
not enter here because only q^qD contribute to (|18| ) 
leading order. We find a color diffuseness parameter c 
fined by j c = Uip c = -D co i or Vp c 

With A ~ g 2 T the analogous result to Eq. (0) is obtain 
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(where we have ignored the minor dependence on N 
Comparing Eqs. (|l6|) and (|2Fj) we see that D co i or 
a s D a = a s Dfi avor . The color-flip mechanism amplif 
the forward collisions so the color cannot diffuse throu 
the quark-gluon plasma as easily as spin or flavor. 

The factor ln(l/a s ) in D co i or has a completely diff< 
ent origin as the one in D flavor or D spin . In D co i or t 
logarithm arises from an integral dq/q over momenti 
transfers from q ~ A ~ g 2 T to q ~ qr> ~ gT as in t 
case of quark and gluon quasiparticle decay rates |J . 
Dfiavor or D sp i„ and the transport coefficient discuss 
in Bpj a similar integral occurs, but with momenti 
transfers from q ~ qn ~ gT to q ~ T. Thus the : 
frared cutoff does not enter these transport coefficien 
In both cases the result is proportional to the logaritl 
of the ratio of the upper and lower limits on the n 
mentum transfer, namely ln(l/(?). The difference in t 
important range of momentum transfers in the two cas 
is due to the absence in the calculation of T) co \ or of t 
extra factor ~ q 2 jT 2 . Therefore small momentum tra] 
fer processes have greater weight in the calculation 
D C oior and quasiparticle relaxation rates than they do : 
Dfiavor, D S pi n and standard transport relaxation rat 
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In the latter cases the factor q 2 /T 2 
reduces the rates by a factor a s . 

Other related transport coefficients are the electrii 
conductivity, a e i, in QED and the corresponding co^ 
conductivity, a co i or , in QCD. Applying a color-elect 
field, E c , to the quark-gluon plasma generates a ca 
current, j c . The color conductivity a co i or = —j c /E c c 
thus be found by solving the Boltzmann equation ana 
gous to the color diffusion process. We find 
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Here D co ior plays the role of the color relaxation tin 
Eq. (21) is the standard result for a plasma except : 



the factor 2/3 which arises because only two-thirds oft 
colors contribute to the currents for a given color fie 
Inserting D co ior from (pfj|) we obtain 
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from quark currents alone. Gluon currents are slower 
due to stronger interactions and will reduce the conduc- 
tivity slightly. This result differs from M by a numerical 
factor only. For comparison the electrical conductivity 
below T<M W is a et ~ T/a\n(l/a) j|. The different 
dependence on coupling constant arises because the ex- 
changed photon does not carry charge whereas the ex- 
changed gluon can carry color. The characteristic relax- 
ation times for conduction are very different in QCD, 
where T color - D color ~ (a s ln(l/a s )T) _1 , as compared 
to QED, where r e ; ~ (a 2 ln(l/cv)T) _1 . Consequently, 
QGP are much poorer color conductors than QED plas- 
mas when T <C Mw for the same coupling constant. 

These surprising results for QCD are qualitatively in 
agreement with those found by Selikhov & Gyulassy [JjJ 
who have considered the diffusion of color in color space. 
They use the fluctuation-dissipation theorem to estimate 
the deviations from equilibrium and find the same color 
non-flip and color-flip terms, which they denote the mo- 
mentum and color diffusion terms respectively, and they 
also find that the latter dominates being infrared diver- 
gent. Inserting the same infrared cut-off they find a color 
diffusion coefficient in color space equal to Eq. (0) 
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Note that this quantity is proportional to the inverse of 
Dcolor as given in Eq. (20). 



The color-flip mechanism is not restricted to QCD but 
has analogues in other non-abelian gauge theories. In the 
very early universe when T ^> M\y ~ 80 GeV, the W r± 
bosons can be neglected and faces the same electroweak 
screening problems as QCD and QED. Since now the ex- 
changed W^ bosons carry charge (unlike the photon, but 
like the colored gluon) , they can easily change the charge 
of, for example, an electron to a neutrino in forward scat- 
terings. Thus the collision term will lack the usual factor 
q 2 as for the quasiparticle damping rate and the color dif- 
fusion. Since SU(2) x U(l) gauge fields should have the 
same infrared problems as SU(3) at the scale of the mag- 
netic mass, ~ e 2 T, we insert this infrared cutoff. Thus we 
find a diffusion parameter for charged electroweak parti- 
cles in the very early universe of order 
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which is a factor a smaller than when T <C Myy- The 
electrical conductivity will be smaller by the same factor 
as well, a e i ~ T/ln(l/a). 

In summary, the flavor, spin and color diffusion coef- 
ficients have been calculated in QCD plasmas to leading 
order in the interaction strength. Color diffusion and 
the quark and gluon quasiparticle decay rates are not 
sufficiently screened and do depend on an infrared cut- 
off of order the magnetic mass, m mag ~ g 2 T; typically 
D co\or ~ a s \n(q D /m mag )T ~ a s ln(l/a s )T. Flavor and 
spin diffusion processes are sufficiently screened by De- 
bye screening for the longitudinal or electric part of the 



interactions and by Landau damping for the transvei 
or magnetic part of the interactions; typically Dj [avor 

D~p in ~ a 2 ln(l/a s )T. As a consequence, color diffusi 
is slow and the QGP is a poor color conductor. In t 
very early universe when T 3> A% exchanges of W p] 
vides charge exchange - a mechanism analogous to co. 
exchange in QCD - and QED plasmas will also be pc 
electrical conductors. 
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